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if) \ ABSTRACT 

■"nI" ■ We report the results of the spatial two-point correlation function £,cc{r) for the new X- 

ray galaxy cluster survey REFLEX, which comprises of 452 X-ray selected clusters (449 
^ ! with redshifts) detected by the ROSAT sateUite during the ROSAT All-Sky-Survey 

^ ■ (RASS). The REFLEX cluster sample is flux limited to 3 x IQ-^^ erg s'^ cm'^ in 

I the ROSAT energy band (0.1 — 2.4 keV) and spans 3 decades in X-ray luminosity 

--^ , (10^^ — 10^^ erg s~^), containing galaxy groups and rich clusters out to a redshift 

• z < 0.3. Covering a contiguous area of 4.24 sr REFLEX is the largest X-ray cluster 

sample to date for which spatial clustering has been analysed. Correlation studies 
Q . using clusters selected on the basis of their X-ray emission are particularly interesting 

$-H ' as they are largely free from the projection biases inherent to optical studies. For 

, the entire flux-limited sample we find that the correlation length (the scale at which 

■ the correlation amplitude passes through unity) tq ~ 20ft-~^Mpc. For example, if a 

^ ' power-law fit is made to ^(r) over the range 4 — 40/i^^ Mpc then rg = 18.8 ± 0.9. 

• i-H , An indication of the robustness of this result comes from the high degree of isotropy 

' seen in the clustering pattern on scales close to the correlation length. On larger 

Jh I scales ^cci'r) deviates from a power-law, crossing zero at ~ 45/i^^ Mpc. From an 

examination of 5 volume-limited cluster sub-samples we find no significant trend of rg 
with limiting X-ray luminosity. A comparison with recent model predictions for the 
clustering properties of X-ray flux-limited samples, indicates that Cold Dark Matter 
models with the matter density = 1 fail to produce sufficient clustering to account 
for the data, while flm — 0.3 models provide an excellent fit. 

Key words: Surveys; Galaxiesiclusters; cosmology: large-scale structure of the Uni- 
verse; X-rays: galaxies 



1 INTRODUCTION 

Clusters of galaxies have been used for many years as tracers 
of the large-scale mass distribution in the universe. As the 
largest gravitationally bound objects their clustering statis- 
tics provide important information on the hierarchical pro- 
cess of galaxy formation enabling estimates to be made of 
the mass fluctuation amplitude and the density paramter 
n (e.g. Mo et al. 1996). The early statistical analyses re- 
lied on the visual cluster compilations of Abell (1958) and 
Abell, Olowin & Corwin (1989). From the redshift surveys 
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of richness-limited subsamples of the Abell catalogue (Bah- 
call & Soniera 1983, Klypin & Kopylov 1983, Postman et 
al. 1992, Peacock & West 1992) it was established that the 
correlation function £,cc{r) followed the form 

M.)^(^)"\ (1) 

on scales < 100/i^^ Mpc with 7 ~ 2 and vq systematically 
5 times higher than the value of ~ 5h~^ Mpc found for 
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galcixies (e.g. Davis & Peebles 1983, Tucker et al. 1997), but 
with a strong dependency on the limiting richness of the 
cluster sample used (see Bahcall 1988). For example, while 
the richness class R > samples give ro ^ 20/i~^Mpc, at 
the higher threshold R > 2 the correlation length rises to 
ro ~ 40/i"^Mpc (Peacock & West 1992). In principle these 
results can be used to place constraints on theoretical mod- 
els of large-scale structure, however the cosmological infor- 
mation they contain is questionable due to the likely ex- 
istence of inhomogencitics (Sutherland 1988, Sutherland & 
Efstathiou 1991) and line-of-sight projection effects (Lucey 
1983, Dekel et al. 1990) artificially enhancing the correlation 
amplitude of Abell-based cluster samples. Strong evidence 
that these effects play a significant role comes from compar- 
ing the amplitude of the correlation function ^(cr, tt) in the 
redshift direction tt of space with the perpendicular direc- 
tion cr. Both rich and poor Abell cluster samples regularly 
fail this isotropy test, showing line-of-sight elongations in the 
contours of ^(cr, tt). These features are consistent with an ar- 
tificial enhancement of the correlation function (Sutherland 
1988, Efstathiou et al. 1992, Peacock & West 1992) although 
physical interpretations have also been suggested (Bahcall et 
al. 1986, Miller et al. 1999). 

The advent of digitised cluster surveys saw a dra- 
matic increase in the homogeneity with which optical clus- 
ter samples could be compiled. Results from both the Edin- 
burgh/Durham Cluster Catalogue (Nichol et al. 1992) and 
the APM survey (Dahon et al. 1992, 1994, Croft et al. 1997) 
demonstrated that for the equivalent Abell richness class 
R>0, clusters found from automated detection algorithms 
have ro ^ 15ft~^Mpc with significantly reduced anisotropies. 
Testing the results of the richer clusters has proved more dif- 
ficult due to the large search volume required to find suit- 
able numbers and the diminishing contrast of distant clus- 
ters against the background of faint galaxies. For example, 
Croft et al. (1997) used 46 APM clusters with richnesses 
equivalent to J? > 2 and found ro — 20 ± 5. 

In recent years attention has focused on cluster samples 
generated on the basis of their X-ray emission. This method 
has enormous advantages for the determination of ^cc over 
the optically compiled cluster catalogues described above: 

• The X-ray emission from a cluster provides a direct 
physical link with the presence of a large gravitational po- 
tential in quasi-equilibrium (e.g. oc M**^^). Thus the sig- 
nature of X-ray emission provides strong evidence that the 
apparent overdensities seen in the optical are gravitationally 
bound structures. 

• The cmissivity of Thermal Bremsstrahlung radiation is 
proportional to the square of the electron number density, 
whereas the optical richness estimates are simply propor- 
tional to the galajcy density. Therefore, at fixed density, the 
contamination in cluster samples resulting from the projec- 
tion of systems along the line-of-sight is intrinsically higher 
in richness-based optical samples compared to X-ray cluster 
catalogues. Furthermore the X-ray emission from clusters 
is concentrated towards the dense central cores which are 
typically ~ 250/i~^ kpc in size - significantly smaller than 
the spatial extent of the galaxy concentration in clusters. 
Both these efi'ects substantially reduce the chance of projec- 
tion effects which, as described above, are thought to plague 
Abell-based samples. 



• The comparatively low internal background of the 
ROSAT Position Sensitive Proportinal Counter (PSPC) and 
the relatively short exposure times in the All-Sky Survey 
means that the X-ray fluxes from clusters at the flux limit 
of REFLEX are photon-noise limited as opposed to back- 
ground limited. This is in contrast to purely optically se- 
lected samples which are forced to have a minimum density 
contrast above a varying background of galaxies before they 
can be detected. 

The first attempts to measure ^cc using X-ray clusters 
were confined to small samples: Lahav et al. (1989) detected 
significant clustering with ro ~ 21h~^ for 7 = 1.8 using an 
all-sky sample of 53 clusters above a flux 1.7 x 10^^^ erg s~^ 
cm"^ (2-10 keV). Nichol, Briel & Henry (1994) used ROSAT 
data for a complete sample of 67 X-ray bright Abell clus- 
ters finding a correlation length ro = 16.1±3.4/i~^ Mpc and 
detecting no significant clustering anisotropy. A more exten- 
sive study using data from the ROSAT satellite carried out 
by Romer et al. (1994) for a nearly complete fiux-limited 
sample of 129 clusters above 1 x 10~^^ erg s~^ cm~^ found 
ro = 12.9±2.2/i"^Mpc, 7 = 1.8±0.4. This study also found 
no evidence of spatial anisotropy in the clustering pattern. 
More recently, there have been two independent estimates 
of ^cc from the 277 X-ray brightest Abell cluster sample 
from the RASS (XBACS, Ebeling et al. 1996). For this sam- 
ple Abadi et al. (1998) suggest ro = 21.1^^ j/i^^Mpc and 
7 = 1.9 from a minimisation procedure using the binned 
correlation data, while Borgani et al. (1999) use a more reli- 
able likelihood analysis finding ro = 26.0't.t'Yh~^ Mpc (here 
error bars are 2a). The anisotropy diagram for XBACS is 
published by Miller et al. (2000) and shows strong Abell- 
typc elongations. 

These X-ray results do not attempt to take account of 
the sky coverage of the parent X-ray survey in the correlation 
analysis. However, X-ray cluster samples generated from the 
RASS (Triimper 1993, Voges et al. 1999) have the advantage 
that the sky coverage is known from accurate information 
on the X-ray flux limit pertaining to amy part of the sky. 

The first attempt to utilise the RASS sky coverage in- 
formation is Moscardini et. al. (2000a), who analyse the spa- 
tial distribution of the clusters in the RASSl Bright Sam- 
ple (De Grandi et al. 1999) using a very simple version of 
the sky coverage based on the first processing of the All- 
Sky Survey. This cluster catalogue is the forerunner to RE- 
FLEX consisting of 130 clusters to a hmit 3 - 4 x 10"^^ 
erg cm~^ defined in the ROSAT hard energy band 
(0.5 — 2.0 keV) and covering an area covering 2.5 sr centred 
on the Southern Galactic Cap. Moscardini et al. (2000a) 
find ro = 21.5+4;!'!"^ and 7 = 2.llgj^ (95.4% errors) with 
a mild dependence of ro on limiting flux and luminosity. The 
REFLEX survey provides the opportunity to substantially 
improve on this result in a number of important respects: 
(i) REFLEX provides more than 3 times the number of X- 
ray clusters over a contiguous area nearly twice as large, (ii) 
All RASS standard analysis source detections are reanal- 
ysed using our own flux determination method, (iii) Due 
account is taken of all exposure variations, in contrast to 
the the RASSl sample which is limited to exposure times 
larger than 150 sees, (iv) The optical identification is done 
in a homogeneous way based on the most comprehensive 
optical data base available for the southern sky. The power 
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Figure 1. The sky coverage for REFLEX as a function of flux 
limit in the ROSAT energy band. The curve is determined from 
the sateflite exposure map, the local hydrogen column density and 
the criterion that at least 10 photons are detected. The dashed 
vertical line shows that REFLEX reaches 3 x 10~^^ erg cm~^ 
for 97.4% of the survey. 

spectrum for REFLEX is presented elsewhere (Schuecker et 
al. 2000), here we concentrate on the corrclatiou function. 

The outhne of the paper is as follows: In Section 2 we 
give a brief description of the REFLEX cluster survey. In 
Section 3 wc discuss the algorithm used to osimatc the the 
correlation function and the results for both the entire RE- 
FLEX catalogue and volume limited sub-samples are pre- 
sented in Section 4. The interpretation of these results in 
terms of structure formation models is discussed in Section 
5. 



2 THE REFLEX SURVEY 

The REFLEX survey represents an objective flux-limited 
catalogue of X-ray clusters in the southern hemisphere south 
of declination +2.5 degs and excluding the region within 
±20 dcg of the Galactic Plane. A further ~ 324 dcg^ of 
sky around the LMC and SMC is removed where X-ray de- 
tection is hampered by the high interstellar absorption and 
crowded star fields. The remaining area covered by the sur- 
vey is 13924 deg or 4.24 sr, representing ~ 34% of the entire 
sky. 

The primary X-ray data for REFLEX originates 
from the second processing of the ROSAT All-Sky-Survey 
(RASS2) using the Standard Analysis Software System 
(SASS) which is based on a maximum likelihood detection 
algorithm. Confirmed RASS2 sources with a likelihood pa- 
rameter of at least 15 and count rate > 0.05 cts s~^ in the 
0.1 — 2.4 keV energy band have already been published in the 
RASS bright source catalogue (Voges et al. 1999). For RE- 
FLEX we use the internal MPE source catalogue totalling 
54076 sources in the study area which allows the inclusion 
of sources with a likelihood > 7. Although some sources will 
be detected at a significance < 3a and not all are real, this 
lower likelihood threshold ensures that the parent catalogue 
is as complete as possible. 

It is well known from previous studies that the RASS 
analysis software is optimised for point-like sources and 



therefore underestimates the flux of extended sources (e.g. 
EbeUng et al. 1996, De Grandi et al. 1997). Therefore we 
have reanalysed all the source fluxes using a growth curve 
analysis method to recover the total flux of extended sources 
with an internal error of between 10 — 20% (Bohringer et al. 
2000a). Note: X-ray count rates are measured in the hard 
band (0.5 — 2.0 keV) then converted to unabsorbed fluxes 
in the ROSAT band (0.1 - 2.4 keV) and the cluster X-ray 
luminosities are determined by an iterative procedure using 
the luminosity-temperature relation of Markcvitch (1998) 
assuming h = 0.5 (Bohringer et al. 2000c in preparation), 
with the values scaled by a factor 0.25 to /i = 1 in this paper. 
Excluding double detections we have 4206 sources above a 
count rate limit of 0.08 cts s^ , which corresponds to a flux 
limit between 1.6 - 2.0 x 10"^^ erg s"^ cm"^. 

The optical identification is based on finding galaxy 
overdensities in concentric rings around X-ray source po- 
sitions using the UK Schmidt J-survey photographic plates 
digitised by COSMOS which reduces the total number of 
cluster candidates to ~ 500 above a flux limit 3 x lO"^'^ erg 
s"^ cm"^ in the ROSAT energy band 0.1 - 2.4 keV. De- 
tails of the optical identiflcation process are given elsewhere 
(Bohringer et al. 2000b). 

To carry out further indentiflcation and obtain red- 
shifts, multi-object (5 — 20 galaxies per cluster) and single- 
slit (2 — 3 galaxies per cluster) spectroscopy was carried 
out on ~ 431 targets as part of an ESO Key Programme 
(Bohringer et al. 1998, Guzzo et al. 1999). This results in 
452 clusters above 3x 10~^^ erg cm~^ in the energy band 
0.1 — 2.4 keV, of which 449 have secure redshifts either from 
our ESO programme or from the literature. About 65% of 
these clusters are in the Abell catalogue while most of the 
others were previously unknown. 

A comprehensive discussion of the contamination and 
completeness statistics in REFLEX is given in Bohringer 
et al. (2000b) and results on the comoving number density 
of clusters in the survey are presented in Schuecker et al. 
(2000). These indicate a completeness well in excess of > 
90% and a contamination by non-cluster X-ray sources of 
less than 9%. We mention a few results here to serve as an 
illustration of the quality of the catalogue: (i) From a search 
for X-ray emission around all AGO and AGO supplimentary 
clusters only 1 cluster with an X-ray flux more than the flux 
limit is not found by the selection process, (ii) Clusters in the 
luminosity range 0.08 — 2. 5x W'^*h~^ erg s"'^ have a constant 
comoving number density of objects and V/Vmax=0.51 ± 
0.01 at the flux limit of the survey. This is consistent with 
the lack of evolution seen in the X-ray luminosity function 
out to at least z ~ 0.3 reported by other surveys {e.g. Burke 
et al. 1997, Ebeling et al. 1997). (ui) Approximately 81% of 
the REFLEX clusters are extended - we searched the RASS2 
database separately for extended X-ray sources flnding only 
a further 8 bona-flde clusters and 5 candidate clusters, 3 of 
which show no obvious optical counterpart and for which 
futher deep imaging is planned. 



3 CALCULATING THE CORRELATION 
FUNCTION 
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3.1 Areal Coverage 



One complication with the RASS is that the sky coverage is 
not homogeneous resulting in about 12% of the REFLEX 
survey region having an exposure time less than half of 
the median exposure time (~ 323 s). This, coupled with 
the varying galactic hydrogen coloumn density, results in a 
variation of the limiting flux of the RASS2 across the sky. 
Although the very low background for the ROSAT PSPC, 
especially in the hard band (0.5 — 2.0 keV), allows the de- 
tection and characterisation of sources with comparitavely 
low source source counts a minimum number is required for 
a safe detection. Fig. |l| shows the resulting effective area 
of the REFLEX survey as a function of flux with the ad- 
ditional criterion imposed of detecting at least 10 photons 
in the hard band. The exposure times of the RASS2 in the 
REFLEX area are sufficient that at a flux limit of 3 x 10~^^ 
erg s"'^ cm~'^ at least 10 photons are detected for 97.4% of 
the REFLEX survey area and hence the number of clusters 
detected with low photon counts is very small - 3.8 clus- 
ters with less than 10 counts are expected in the survey and 
only 1 is detected. For a more conservative requirement of 
at least 30 photons for each source the sky coverage falls to 
78% (see Bohringer et al. 2000b). 



3.2 Correlation Estimator 

In all computations we use the estimator 



l + ^cc(r) =4 



{DD){RR) 
(Di?)2 ' 



(2) 



where DD stands for the number of distinct pairs in 
the data, RR stands for the number of distinct pairs in the 
random catalogue and DR represents the number of cross 
pairs. The factor of 4 in this expression accounts for the fact 
that while the number of distinct pairs in a large catalogue 
of size n, say, is ~ /2, the number of cross-pairs between 
two different catalogues each with n entries is ~ and these 
are all distinct. Since the correlation function is defined in 
terms of the total number of pairs, the number of distinct 
DD and RR pairs must each be multiplied by 2 to obtain 
the total numbers, hence the factor of 4. This estimator has 
been shown by Hamilton (1993) to be the most robust for 
datasets which may be sensitive to the chance location of 
strong clustering close to the sample boundary. 

In principle the estimator in equation ^ can be gener- 
alised to include an arbitrary weighting function. For calcu- 
lating the correlation function of galaxies from magnitude- 
limited samples, the variance in the estimate of ^(r) on large 
scales is minimised if w{ri,T) is 



1 + 4TTnDJ3{r)(j}{ri 



(3) 



where is the distance of an object from the origin, 
T is the distance separating two objects, 4){r) is the survey 
selection function, no the mean space density of objects 
and Jiir) — drr^S,{r) (see Saunders, Rowan- Robinson 
& Lawrence 1992, Fisher et al. 1994, Guzzo et al. 2000). 
The physical interpretation of the term A-KnoJiij) in the 
weighting scheme is that it represents something like 'the 
mean number of objects per clump'. For galaxies this num- 
ber is large on small scales giving equal volume weighting to 




Figure 2. Aitoff plot of 449 REFLEX Clusters with redshift. 

the pairs (w oc l/(j)(r)). By contrast, for clusters the weight- 
ing term AimoJiij) is always small compared to unity on 
scales of interest and consequently we assign equal weight to 
all pairs in the calculation of the cluster correlation function. 

We calculate spatial separations using the formula for 
comoving coordinate distance (ri); 



c 



(go^) + ((l-go)(l-((2go2) + l)^/')) 



(4) 



adopting the cosmology Ho — IQQh^^ Mpc, Q,m — 1.0 
& S^A = 0.0, along with the cosine rule to determine angular 
separations. 



3.3 Random catalogues 

The random catalogues are constructed over the REFLEX 
survey area using a Monte-Carlo technique which incorpo- 
rates knowledge of the flux limit in cells of size ~ 1 square 
degree. To begin with it is assumed that the observed num- 
ber count distribution of X-ray clusters, LogN-LogS, is well 
fitted by a simple lower-law: 



7V(> S) = AS- 



(5) 



For the purposes here we adopt the value a — 1.35, 
consistent with the REFLEX number counts (see Bohringer 
et al. 2000b) and those of RASSl, the precursor survey of 
REFLEX (De Grandi et al. 1999). Small changes to the value 
of a does not alter the outcome of the results. If we assume 
that the accumulative distribution P defined as 



iV(< S) 



Ntotal{> Slim) 



= P 



(6) 



is uniformly distributed in the range — > 1, then the 
distribution of cluster X-ray fiuxes 5* selected at random 
above Sum is given by 



Sumii^py 



(7) 



We select a cluster at random within the allowed bor- 
ders of the REFLEX survey and then use eqn. |^ to assign 
it a fiux. We then test whether the cluster falls above or 
below the fiux limit for that region of the REFLEX survey 
based on the local values of exposure time and the inter- 
stellar hydrogen column density (Dickey & Lockman 1990, 



© 0000 RAS, MNRAS 000, 000-000 



The spatial correlation function of REFLEX clusters 5 




Figure 3. Aitoff plot of 1000 random points generated with the 
REFLEX sensitivity map (flux limit 3 X 10~^^ erg cm~-^ with 
a minimum of 10 photon counts) and mask. 
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Figure 5. Histogram of the 449 REFLEX clusters with redshifts. 
The solid line is the estimated density distribution using a Gaus- 
sian kernel of width 5600 km . 
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Figure 6. Histogram of the 399 REFLEX clusters with > 
3 X lO** erg s~^. The dashed curve is the estimated redshift dis- 
tribution of this sample using a Gaussian kernel of width 5600 
km s~^, as for Fig. |^. The solid curve is the redshift distribution 
estimated by integrating the X-ray cluster luminosity function. 



Figure 4. Histogram of REFLEX clusters as a function of galac- 
tic longitude and latitude. The solid line is the prediction from a 
random distribution of points convolved with the REFLEX sen- 
sitivity map (flux limit 3 X lO"^'^ erg s~^ cm"'^ with a minimum 
of 10 photon counts) and mask. 



Stark et al. 1992). We set Su^^ = 3 x 10"^^ erg s"^ cm-^ 
which is the cut-off flux limit for the entire survey. An ex- 
ample of a random catalogue with 1000 points generated in 
this way is shown in Fig. ^. To demonstrate reliability of 
the random catalogues the histogram in Fig. ^ shows the 
number of REFLEX clusters as a function of Galactic longi- 
tude and latitude compared to that of a random catalogue 
generated using the REFLEX survey sensitivity map. The 
random catalogues used in the determination of ^cc contain 
typically 100,000 points. 

We use two methods to assign each random point a 
redshift. (i) Redshifts are drawn from the distribution of 
REFLEX clusters smoothed with a Gaussian kernel. This 
method allows the redshift selection function to be estimated 



without prior knowledge of the underlying density distribu- 
tion of clusters and is used in almost all previous determi- 
nations of the cluster correlation function. For REFLEX we 
use a Gaussian of width 5600 km s~^ - the optimum value 
depends on the space density of clusters and is constrained 
by the need to follow the redshift distribution accurately 
enough while not removing large-scale clustering. The ex- 
act figure used is generally not critical, with values in the 
literature ranging between 4000 — 8000 km s~^. (ii) The sec- 
ond method, which we apply to luminosity limited samples, 
uses the X-ray cluster luminosity function to generate the 
expected number of clusters at each redshift. Assuming a 
Schechter function of the form 



n{L)dL = Aexp{L/L^,){L/L^ 



(8) 



where n(L) is the number density of clusters per lu- 
minosity interval, then for particular values of a and i*, 
we can integrate n{L)dL above Lum to determine the num- 
ber density of clusters at each redshift The value of 
Liim at each redshift is found from the flux limit (fixed at 
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3 X lO^^^erg cm~^). The expected number of clusters in 
each redshift interval dz is then simply rj(z) * dV[z), where 
dV{z) is the comoving volume element at redshift z. Each 
random point generated in the area of the survey is thus 
assigned a random redshift weighted by the expected num- 
ber of clusters based on eqn. ^ We have used the values of 
a = 1.61, L* = 6.04 x lO"" (erg s"^), A = 3.04 x 10"** (10'*'' 
erg s~^)~'^ and Ho = 50 km s~^ Mpc~^, which are appro- 
priate to the REFLEX sample (Bohringer et al. 2000c, in 
preparation), although using our previous luminosity func- 
tion parameters from De Grandi et al. (1999) gives identical 
results. The redshift distribution for the REFLEX sample 
along with the smoothed version using method (i) is shown 
in Fig. 1^, while Fig. ^ shows a comparison between the red- 
shift distribution of both methods for the luminosity sub- 
sample of 399 clusters limited to L^, > 3 x lO^^erg s~^. 
We prefer method (ii) for the luminosity sub-samples as it 
makes no prior assumptions regarding the scale of the clus- 
tering and avoids the need to smooth the data, however in 
practice we found no significant difference in the correlation 
functions resulting from the two methods. 




r (h ' Mpc) 

Figure 7. The correlation function for the REFLEX survey of 449 
clusters. The error bars on each point are derived from bootstrap 
statistics. The solid line shows the result of the likelihood analysis 
(ro = 18.8, 7 = 1.83) fitting a power law over the range 4 — 40/i~"'^ 
Mpc. 



3.4 Maximum Likelihood Determination of ro and 

7 

In calculating the best-fit power-law for the correlation func- 
tion from samples of ~ 100 clusters there has traditionally 
been one of two methods adopted. The first is to calculate er- 
rors for ^cc based on estimates from pair counts binned into 
~ 10 coarse intervals, which are usually spaced logarithmi- 
cally out to ~ lOO/i"^ Mpc {e.g. Bahcall et al. 1983, Dalton 
et al. 1992, Nichol et al. 1992). On the grounds that each 
bin contains a large numbers of pairs, the best-fit power-law 
is estimated using the statistic. The danger with such an 
approach is that the resulting estimates of parameters de- 
scribing the power-law (ro, 7) will then depend on the precise 
details of the binning. In order to overcome this limitation 
we adopt the second of the two methods referred to above 
and maximise the likelihood L that the model correlation 
function produces the measured number of cluster pairs at 
a given separation (Croft et al. 1997, Borgani et al. 1999, 
Moscardini et al. 2000a). The likelihood estimate is based 
on Poisson probabilities, such that 

JV 

L = l[e-^^''/u\, (9) 

i=l 

where u is the observed number of cluster-cluster pairs 
in a small interval dr and fj, is the expected number in the 
same interval calculated using Hamilton's estimator (eqn.|^. 
As long as the number of random points is kept large enough 
to avoid a zero in the denominator of eqn. ^ dr can be 
made arbitrarily small, so as to ensure the final results are 
independent of the bin size. In practice we used ~ 7000 
bins between 5 — 100 Mpc, which resulted in either a or 1 
cluster-cluster pair in almost all bins. 

The associated errors on the correlation function are 
usually calculated from 'Poisson' statistics. In the case of 
large bin intervals errors are computed from the formula 

(1 +Ccc(r)) 

S^cc[r) = -= — , (10) 



where Ncc is the number of distinct cluster pairs in the 
bin centred at separation r. In the case of a maximum likeli- 
hood determination, such as that used here, confidence lev- 
els can be defined as S{rt,est,'ybest) — 'S'(ro, 7), where S is the 
usual S = —2lnL, assuming that AS is distributed like x^- 

Both these methods are likely to produce underesti- 
mates of the true dispersion as the use of Poisson statistics 
assumes that the pair counts are independent of each other, 
which is clearly not the case. An estimate of the true dis- 
persion in the correlation function, which tries to account 
for cosmic variance, can be made either by applying a boot- 
strap resampling of the real data {e.g. Ling, Frenk & Barrow 
1986; Mo, Jing & Borner 1992) or by carrying out numerical 
simulations based on plausible cosmological models (Croft 
& Efstathiou 1994, Croft et. al. 1997). Both methods pro- 
duce similar results indicating that the real errors are prob- 
ably 1 — 2 times larger than the Poisson-based estimates. 
We confirmed this result for our sample by generating mock 
catalogues from bootstrap resampling of the data and cal- 
culating the best-fit ro and 7 values using the likelihood 
method described above. The ratio of the error for ro from 
the variance between the bootstrap samples and the Poisson 
error is between ~ 1.5 — 2.0 for all luminosity subsamples. 
Unless stated otherwise, in the results which follow we quote 
the la likelihood errors on the values of ro and 7. 



4 RESULTS 

The 5cc(r) for the REFLEX survey of 449 clusters is shown 
in Fig. |. Afit was made to the correlation function assuming 
a single power law over the range 4 — 40h~^ Mpc using the 
likelihood analysis described in Section 3.4. Fig. ^ shows the 
corresponding joint constraints resulting from this analysis. 
The best-fit value for the power-law parameters are ro — 
18.8 ±0.9 and7 = 1.83+o:o8. If points are included on larger 
scales then the slope steepens, e.g. fitting between 4— 100/i~^ 
Mpc gives 7 ~ 2.35 and ro — 16.25. The inability of a single 
power law to adequately describe the correlation function is 
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Figure 8. The probability contours {la, 2a, Str) for the best fit 
ro and 7 from the likeUhood analysis of the REFLEX survey over 
the range 4 - 40/i~l Mpc. 



Figure 10. The correlation function amplitude ro plotted against 
limiting X-ray luminosity defined in the ROSAT energy band 
(0.1 — 2.4 keV). The errors correspond to la from the likelihood 
analysis. 
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Figure 9. The X-ray luminosity (0.1 — 2.4 keV) vs redshift for 
the REFLEX sample of 449 clusters. Also shown are the 5 vol- 
ume complete sub-samples used to examine the variation of the 
correlation length with limiting X-ray luminosity {Lu^. 
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Figure 11. The correlation function amplitude in the range — 
20h~^ Mpc plotted against limiting X-ray luminosity defined in 
the ROSAT energy band (0.1 - 2.4 keV). The error bars are la 
based on eqn. 



further reflected in the zero crossing of ^cc at 4:5h~^ Mpc 
(see Section 7.1). 

In order to investigate the dependency of ro with 
X-ray luminosity we also calculated ^cc{r) for 5 volume- 
limited X-ray sub-samples with luminosity thresholds 
0.08,0.18,0.3,0.5,1.0 in units of lO^^h-^ erg s'V Fig. | 
shows the distribution of luminosity with redshift for the 
REFLEX sample along with the regions corresponding to 
the 5 subsamples. Due to the significant covariance be- 
tween ro and 7 this investigation has been carried out 
with 7 fixed at 2.0. The correlation results for the sub- 
samples are presented in Table ^ and Fig. These in- 
dicate no significant positive trend of ro vs Lurn, with the 
highest measured ro occuring at intermediate luminosities 
(iiim > 0.3 X lO**/!"^ Mpc). Beyond this point the sta- 
tistical errors increase rapidly. To test the reliability of the 
parameter ro as an indicator of how clustering changes with 
X-ray luminosity, we calculated the average correlation func- 
tion amplitude for the 5 volume-limited sub-samples over the 
range of separations — 20h~^ Mpc. The result, shown in 
Fig. is in good agreement with the trend of ro vs Liim- 




Figure 12. Contours of constant ^(a, n) for the REFLEX survey. 
The contour values used are 4.0, 2.0, 1.0 and then 0.8 to -0.4 in 
steps of 0.2 
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Table 1. Estimates of ro as a function of limiting X-ray lumi- 
nosity for volume limited cluster subsamples. Values of ro are 
calculated at 7 = 2.0 



L,,„{1044h-2) ergs 



Number 



1.0 

0.5 

0.3 

0.18 

0.08 



67 
101 
108 
84 
39 



22.9 

25.8 

31.1 

25 

24.8 



+7.3 
7.7 
+3.2 
3.3 
+2.0 
2.1 
1.9 
2.0 
+2.5 
2.5 



We have examined the isotropy of the clustering sig- 
nal for the REFLEX survey by plotting contours of ^(cr, tt), 
where it = \ri — r2\ is the line-of-sight separation, with ri,2 
determined from eqn. ^, and a = (s^ — tt^)^'''^ is the perpen- 
dicular component of the cluster separation s. As discussed 
in the introduction, elongations of the contours in the red- 
shift direction compared to the perpendicular direction for 
scales ~ 2Qh~^) Mpc are a feature of some optical cluster 
catalogues - typically with a ratio ~ 4 : 1 for redshift sam- 
ples based on the Abell catalogue (e.g. Postman et al. 1992). 
Fig. represents the corresponding plot for the REFLEX 
clusters and indicates that unlike optical surveys, the f (a, n) 
contours are close to being completely concentric on scales 
close to the correlation length. 



5 DISCUSSION 

The determination ^cc from the REFLEX survey can be 
compared with similar determinations for other X-ray clus- 
ter samples. Our results of ro = 18.8 and little dependency 
of ro on X-ray luminosity are broadly consistent with the re- 
suhs of XBACS (Borgani et al. 1999) and RASSl (Moscar- 
dini et al. 2000a). The result presented by Romer et al. 
(1994), hereafter R94, for a sample of 128 clusters above 
1.0 X 10"^^ erg s"^ cm' 



in a 3100 deg'^ area centered on 
the SGP, suggests a correlation length ro = 13 — 15/i~^ 
Mpc, smaller than any of the other determinations from X- 
ray samples. In addition to the fainter flux limit, the R94 
study differs from REFLEX in 2 further ways which in prin- 
ciple could affect the result: (i) the cluster sample was based 
on a reduction of the all-sky-survey using the ROSAT Stan- 
dard Analysis Software, which has subsequently been re- 
vised (ii) the correlation analysis performed by R94 did not 
include the sample sky coverage corresponding to the SGP 
region under study. We have investigated the origin of a 
possible systematic difference between R94 and REFLEX 
by repeating our correlation analysis on the 109 REFLEX 
clusters lying within the R94 SGP area of sky, defined by the 
boundaries 22hr < RA < 3hr, -50° < dec < 2°,|6| > 40°. 
The resulting power-law fit to the correlation function out to 
< lOO/i"^ Mpc gives ro = I2.'dt\l 7 = 2-Olo:t smaller than 
the REFLEX amplitude of 18.8 ± 0.9 and very close to the 
original SGP result of ro = 12.9±2.2/i"^ Mpc, 7 = 1.8±0.4 
found by R94 fitting over the same range. This suggests that 
the difference between the REFLEX and the SGP result is 
most likely due to the superior statistical sampling of RE- 
FLEX which represents a 4-fold increase in survey area over 
the SGP region while probing to a similar redshift. 

Miller et al. (2000) analyse the ^{a,n) diagram for 




Figure 13. The REFLEX correlation function compared with a 
range of Cold Dark Matter models for which Om = 1 and Ha = 0, 
using the same flux limit and sky coverage as REFLEX - taken 
from Moscardini et al. (2000b). Error bars in this figure are la 
and based on bootstrap resampling. 



a number of X-ray cluster samples. In their analysis 
the XBACS clusters show very strong elongations around 
^(ct, tt) ~ 1 in the redshift direction and a similar anisotropy 
is present in other X-ray confirmed Abell cluster samples. 
The RASSl sample shows a much weaker anisotropy over 
the same scale. On the basis of this Miller et al. (2000) argue 
that clustering anisotropy is a ubiquitous feature of X-ray 
cluster samples which demonstrates that the anisotropics 
are real. However, the absence of any significant anisotropy 
in the contours of ^((t, tt) for the REFLEX survey shown 
in fig. |l^ indicates that this is not the case. This is the 
strongest indication yet that elongations seen in other cata- 
logues are spurious and justifies the claim, first pointed out 
by R94, that X-ray selected surveys do not suffer from sig- 
nificant projection biases. As with optical studies based on 
the Abell catalogue, the presence of strong elongations close 
to the scale of ro in the XBACS bring the accuracy of the 
clustering signal derived from this sample into question. 



5.1 Comparison with Cosmological Models 

Predictions for the clustering properties of X-ray selected 
clusters from a number of surveys, including REFLEX, have 
recently been made by Moscardini et al. (2000b). In these 
predictions the structures on a given scale are assumed to 
evolve by hierarchical merging of smaller units and instan- 
taneous merging on cluster scales. The comoving mass func- 
tion of haloes is computed using the Press-Schechter (1974) 
technique but incorporating more recent corrections which 
improve the comparison of Press-Schechter with numerical 
simulations (e.g. Sheth & Tormen 1999). The link between 
X-ray luminosity and mass of the hosting dark matter halo 
begins with the empirical relation between gas temperature 
T and X-ray luminosity Li,oi 



(11) 



with A = 4.2 & 13 = 1/3, which is a good approxi- 
mation for clusters (e.g. David et al. 1993; White, Jones & 
Forman 1997; Markevitch 1998). The parameter v describ- 
ing the evolution of the T — Lboi relation is constrained by 
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Figure 14. The REFLEX correlation function as for Fig. 
but compared with two open CDM models for which Q,m = 0.3, 
Qa = 0.0 (OCDM) and = 0.3, Ha = 0.7 (ACDM), also from 
Moscardini et al. (2000b). 
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Figure 16. The REFLEX correlation function at low amplitude. 
The error bars are calculated from eqn. hol Also shown is the 
best power-law fit over the range 4 — 40/i Mpc corresponding 
to ro = 18.8 and 7 = 1.83. Predictions for the open and A- 
dominated CDM models from Moscardini et al. (2000b) are also 
presented. 
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Figure 15. Comparison of the REFLEX correlation length tq 
from the likelihood analysis with predictions as a function of lim- 
iting X-ray flux for a range of CDM-type cosmological models 
(see text for details). The REFLEX point is plotted at 3 X 10" 
erg s~^ cm~^. The likelihood error bar of 0.9 on ro has been 
increased by a factor 1.5 to 1.35 here to account for the extra 
contribution due to cosmic variance as described in Section 3.4. 



Moscardini et al. (2000b) using the X-ray cluster number 
counts over the range 5 x 10"^'^ — 3 x 10~^^ erg s~^ cm~^ 
(0.5 - 2.0 keV) taken from the RASSl Bright Sample (De 
Grandi et al. 1999) and the fainter ROSAT Deep Cluster 
Survey Rosati et al. (1998). It is possible to convert the 
temperature estimates from eqn. fOI to halo mass assuming 
a virial isothermal gas distribution and spherical collapse 
(e.g. Eke, Cole & Frenk 1996) . Finally, in order to make pre- 
dictions for the correlation function of the REFLEX survey, 
Moscardini et al. (2000b) incorporate the actual sky cover- 
age of the survey shown in Fig. hi for the passband 0.1 — 2.4 
keV in an identical manner to the procedure used for calcu- 
lating £,cc{r) using the REFLEX data described in Section 
3.1 above. 

The behaviour of the cluster correlation function for 
a range of popular cosmological models based around cold 
dark matter (CDM) are shown in Fig. ^ and Fig. |l|. The 



model predictions are taken directly from Moscardini et al. 
(2000b) and represent: Standard (SCDM), r (rCDM) and 
tilted (TCDM) models, all with fi™ = 1 and ^a = 0.0; along 
with an open model Q.m = 0.3 (OCDM) and a A-dominated 
model with Sim = 0.3 and SIa = 0.7. An examination of 
these figures reveals clear evidence for an inconsistency be- 
tween all Q.m = 1 models and the REFLEX correlation func- 
tion, with a significantly better fit to the data for open or A- 
dominated cosmologies. This is further illustrated in Fig. ^ 
which shows the Moscardini et al. (2000b) predictions of cor- 
relation length with limiting X-ray fiux. While the Einstein- 
de Sitter models predict ro = 11 — 13/i~^Mpc, both the 
OCDM and ACDM models predict ro ~ 20/i"^Mpc. These 
results are in agreement with the analysis of the X-ray clus- 
ters based on the RASSl Bright Sample (Moscardini et al. 
2000a) and the digitised optical surveys (e.g. Croft et al. 
1997). 

Confirmation of the general conclusions on the form of 
the cosmological power spectrum comes from the behaviour 
of of ^cc on large scales. In Fig. |l6| we show the REFLEX 
^cc(r) at low amplitude which shows a positive clustering sig- 
nal out to at least 40/i"^ Mpc (^(30-40) = 0.32±0.08), with 
a zero-crossing ~ 45/i"^ (^(40 - 50) = -4.7 x 10"^ ± 0.05). 
On larger scales the amplitude remains slightly negative 
(C(50- 100) = -0.07 ±0.02). Also shown is the curve repre- 
senting the power law with ro = 18.8/i~^ Mpc and 7 = 1.83, 
along with predictions from Moscardini et al. (2000b) for the 
OCDM and ACDM models. Since the SCDM model predicts 
a zero-crossing near r ~ 33/i~^Mpc (Klypin & Rhee 1994), 
the REFLEX data again support the findings of other clus- 
ter surveys that models with more power than SCDM are 
required to adequately fit the large-scale Ccc(r). Generally 
for CDM-like models (with n = 1 for the primordial spec- 
tral index) the first zero-point occurs at r ~ 16.5(^2™,^^)"^ 
for a vanishing baryon fraction (see Klypin & Rhee 1994). 
For the particular parameterisation used by Moscardini et 
al. (2000b) the OCDM and ACDM models remain positive 
until ~ 80/i~^ Mpc, however the small clustering amplitude 
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in the data between ~ 45 — lOO/i ^ Mpc seen in Fig. 
prohibits any definitive comparison of the zero crossings. 



6 SUMMARY 

Cataiogues of galaxy clusters based on their X-ray emission 
provide a powerful tool for studies of large-scale structure. 



We present tne spatial correlation tunction ot tne itJii^ L/£jA 



cluster kurvey, whldi cuiihIhlh 01449 X-ray ymluliig t;liiyLei'H 
above a flux limit 3 x 10^'' erg s~^ cm~^ and covering a 
contiguous area of 4.24 sr in the southern hemisphere. The 
advantages of X-ray selection, combined with the increased 
statistics and high completeness of REFLEX enable a signifi- 
cant step to be taken in establishing the clustering properties 
of clusters in the local universe. Over the scale 4 — 40/i~^ 
Mpc we find a correlation amplitude ro = 18.8 ± 0.9 and 



power law index 7 = 1.83_ 



for the entire survey. The 



high degree of isotropy in the correlation function demon- 
strates that systematic projection effects are not present in 
the data. By analysing volume-limited sub-samples we find 
no significant trend of clustering amplitude with X-ray lumi- 
nosity. Comparing the REFLEX ^cc results with predictions 
from various CDM-type models which incorporate directly 
the areal coverage of REFLEX, Q.m — 0.3 models provide 
an excellent fit, while fim = 1 & = models fail to 
provide enough large-scale power. Finally, it is intriguing to 
note the concensus emerging between clustering studies and 
the lack of evolution in the abundance of X-ray clusters (e.g. 
Burke et al. 1997, Collins et al. 1997, Borgani et al. 1999, 
Henry 1997, Nichol et al. 1999), which also indicates that 
the Einstein de-Sitter universe is in trouble. 
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